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Abstract 

We study the impact of neutrino masses and mixings on LFV processes within the 
context of the supersymmetric seesaw scenario, where the CMSSM is extended by three 
right-handed (s)neutrinos. A hierarchical spectrum is considered for both heavy and 
light neutrinos. We systematically analyse the interesting relation between the leptonic 
mixing angle #13 and LFV muon and tau decays, namely lj — * 1%^ and lj — > 3U, and 
discuss the interplay with the other relevant parameters. We require compatibility 
with low energy neutrino data, bounds on both LFV decays and charged lepton elec- 
tric dipole moments, and impose a successful baryogenesis via thermal leptogenesis. 
Particular emphasis is given to the implications that a future #13 measurement can 
have on our knowledge of the heavy neutrino sector. 
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1 Introduction 



The impressive experimental evidence of neutrino masses [1] has lead to the first clear signal 
of physics beyond the standard model (SM). One of the simplest extensions of the SM that 
allows to naturally explain the smallness of the neutrino masses (without excessively tiny 
Yukawa couplings) consists in incorporating right-handed Majorana neutrinos, and imposing 
a seesaw mechanism for the neutrino mass generation [2, 3]. The seesaw mechanism offers in 
addition the interesting possibility of baryogenesis via leptogenesis [4] . Within the framework 
of leptogenesis, the observed baryon asymmetry of the Universe (BAU) is explained by the 
out-of-equilibrium decays of the same heavy right-handed neutrinos which are responsible for 
the suppression of the light neutrino masses. The scale of new physics is naturally introduced 
by the heavy right-handed neutrino masses which, for the simplest case of just one right- 
handed neutrino, and assuming neutrino Yukawa couplings, Y u , of 0(1), typically lies close 
to 10 14 GeV. 

Supersymmetric (SUSY) extensions of the SM, including three right-handed neutrino 
superfields, are well motivated models which can accommodate a seesaw mechanism, and at 
the same time stabilise the hierarchy between the scale of new physics and the electroweak 
(EW) scale. One of the most striking phenomenological implications of SUSY seesaw models 
is the prediction of sizable rates for lepton flavour violating (LFV) processes [5]. Assuming 
Y v ~ 0(1) and that the scale of soft-SUSY breaking is of the order (or below) 1 TeV, the 
radiatively induced SUSY corrections driven by the neutrino Yukawa couplings lead to rates 
for the LFV observables which are many orders of magnitude larger than those expected 
from the SM seesaw. Even though this holds irrespective of the chosen mass pattern for the 
right-handed neutrinos, it has been shown that when compared to the degenerate case [6, 7] 
the hierarchical scenario leads to larger LFV rates, which may even be within the reach of 
current experimental bounds [8-20]. In this sense, the lj — > k'-f and lj — * 3Zj (i ^ j) lepton 
decay channels, as well as /i — e conversion in heavy nuclei, are among the most interesting 
processes [21]. Experimentally, the most promising decay is the /i — > process, which 
exhibits the most stringent present bounds, and offers a significant improvement regarding 
the future sensitivity. Furthermore, in the presence of complex neutrino Yukawa couplings, 
one can also construct from the latter LFV decays interesting observables, which are sensitive 
to CP violation in the neutrino sector. For instance, one can build T- and P-odd asymmetries 
in \i — > e7 and /i — > 3e decays, which were addressed in [22,23]. 

In addition to the large number of parameters of the minimal supersymmetric standard 
model (MSSM), the seesaw mechanism introduces 18 new parameters in the neutrino sector. 
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As a first step to simplify the analysis of the LFV rates in a SUSY seesaw model, we choose 
to work in the so-called constrained MSSM (CMSSM), assuming universality of the soft- 
SUSY breaking parameters at the scale of gauge coupling unification, M X - This allows 
to reduce the unknown parameters in the SUSY sector to the five usual parameters of a 
minimal supergravity (mSUGRA) framework. Moreover, regarding the neutrino sector, we 
will assume a hierarchical spectrum for both light and heavy neutrinos. 

Among the various seesaw parameters, those connected to the light neutrino sector can 
in principle be derived from low-energy neutrino data, while those associated with the heavy 
neutrino sector, in particular the values of the masses, are a priori clearly unreachable. How- 
ever, and given the fact that both light and heavy neutrinos enter in the determination of the 
LFV rates (via the Yukawa interactions), a powerful link between the low- and high-energy 
neutrino parameters can be obtained from these LFV processes. From the requirement of 
compatibility with both current LFV bounds and low-energy neutrino data, one can then 
extract information on the heavy neutrino sector, thus providing an indirect access to the 
heavy neutrino parameters. In the presence of additional CP phases (other than those asso- 
ciated to the light sector), and assuming that BAU is generated from thermal leptogenesis, 
one can obtain a further insight on the heavy neutrino parameters. More specifically, one 
can obtain a lower bound on the mass of the lightest right-handed neutrino, which in turn 
translates into lower mass bounds for the other heavy states. 

Here we address the subject of how to extract information on the unknown SUSY-seesaw 
parameters from both the analysis of LFV decays and the requirement of successful BAU. 
As already said, we restrict ourselves to the scenario of hierarchical heavy neutrinos, which 
leads to the most interesting predictions. There are studies also addressing the same subject, 
within similar SUSY seesaw frameworks and hierarchical neutrino scenarios. In particular, 
some constraints on the heavy neutrino and SUSY sectors have been considered by [13-20] 
and further implications regarding low-energy CP violating observables were studied in [22- 
25]. In addition, it has been noticed that not only are the LFV branching ratios (BRs) 
sensitive to the heavy neutrino parameters, but there is also a potentially relevant role being 
played by the yet undetermined low-energy neutrino parameters. The latter include the 
overall light neutrino mass scale and CP violating phases. Concretely, the effects of this 
scale and of the two Majorana phases on the radiative LFV decays have been analysed 
in [26]. 

An even more intriguing situation occurs regarding the sensitivity of the LFV rates to 
$13, which is our main interest here. Although it is not expected to be of relevance in the 
context of the LFV phenomenology, the few studies regarding the dependence of LFV rates 



3 



on #13 [15, 18] have revealed that, for some specific seesaw cases, the LFV muon decays 
indeed exhibit a strong sensitivity to this parameter. In particular, the dependence of the 
BR(/i — > e 7) on #13 was noticed [15] in the context of SUSY grand unified theories (SUSY- 
GUTs), and in [18] within the CMSSM. In the latter, it was further emphasised that in 
addition to /i — > cy other LFV decays as /i — > 3e are also very sensitive to 9 13 . Both these 
studies assumed a simple scenario where no additional mixing, other than that induced from 
the low-energy mixing angles %, was present in the neutrino sector. 

In this work we systematically explore the sensitivity of LFV processes to 9 1S in a broader 
class of SUSY seesaw scenarios, with different possibilities for the mixing in the neutrino 
sector, and we incorporate in our analysis the requirement of BAU via thermal leptogenesis. 
Specifically, we conduct a comprehensive and comparative study of the dependence on 9i 3 
in all the following decay channels: fj,~ — > e~ 7, /i - — > e + , t~ — > e~ 7, r~ — > e~ e~ e + , 

t~ — > /i _ 7 and t~ — > fi~ fi~ We will show here that various of these channels indeed 
offer interesting expectations regarding the sensitivity to #13. In our analysis, we work 
in the context of a CMSSM extended by three right-handed neutrinos and their SUSY 
partners, and use the requirement of generating a successful BAU in order to constrain 
the explored seesaw parameter space. Our main motivation to perform the present study 
has been triggered by the potential measurement of #13, as suggested by the experimental 
program of MINOS and OPERA, which claim a future experimental sensitive of #13 < 8° [27] 
and #13 < 7° [28, 29], respectively. With the advent of other experiments, like Double Chooz 
and T2K, the expected sensitivity will be further improved to #13 < 4° [30] and #13 < 2° [31]. 
An ambitious program to push the sensitivity to less than 1° is envisioned by Neutrino 
Factory [32,33] and/or Beta Beam [34,35] facilities. 

Our ultimate goal is to explore the impact of a potential #13 measurement on the LFV 
branching ratios, which together with the current and future experimental bounds (measure- 
ments) on the latter ratios, could lead to a better knowledge (determination) of the heavy 
neutrino parameters. 

Our work is organised as follows. In Section 2, we present the SUSY seesaw scenario, de- 
scribing the seesaw mechanism for the neutrino mass generation, and discussing how flavour 
mixing in the slepton and sneutrino sectors arises in this context. We further address the con- 
straints on the seesaw parameters from the requirement of generating a successful BAU via 
thermal leptogenesis, and from imposing compatibility with experimental data on charged 
lepton electric dipole moments (EDMs). In Section 3, we explore in detail how the several 
parameters affect the theoretical predictions for the LFV rates, and whether the former can 
modify the sensitivity of a given LFV observable to 9 13 . Section 4 is devoted to the discus- 
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sion of the hints on SUSY and seesaw parameters which can be derived from a potential 
measurement of 6> 13 and LFV branching ratios. Finally, our conclusions are summarised in 
Section 5. 



2 LFV within the SUSY seesaw model 

In what follows, we first present the SUSY seesaw scenario within the CMSSM, then proceed 
to describe how LFV processes arise in this framework and finally discuss the implications 
regarding BAU and charged lepton EDMs. 

2.1 The SUSY seesaw scenario 

The leptonic superpotential containing the relevant terms to describe a type-I SUSY seesaw 
is given by 

W = N C Y V LH 2 + PYiLH! + \N c m M N\ (1) 

where N c is the additional superfield that contains the three right-handed neutrinos v R . and 
their scalar partners vr v The lepton Yukawa couplings Y t;l/ and the Majorana mass m M are 
3x3 matrices in lepton flavour space. From now on, we will assume that we are in a basis 
where Yi and Tom are diagonal. 

After EW symmetry breaking, the charged lepton and Dirac neutrino mass matrices can 
be written as 

m { =Y l v x , m D = Y v v 2 , (2) 

where Vi are the vacuum expectation values (VEVs) of the neutral Higgs scalars, with v 1 ( 2 ) = 
v cos(sin)/5 and v = 174 GeV. 

The 6x6 neutrino mass matrix is given by 

iir = ( ml V (3) 

y m D m M J 

The eigenvalues of M v are the masses of the six physical Majorana neutrinos. In the seesaw 
limit, the three right-handed masses are much heavier than the EW scale, ^ v, and 
one obtains three light and three heavy states, ui and N iy respectively. 

Block-diagonalisation of the neutrino mass matrix of Eq. (3), leads (at lowest order in the 
(m D /m M ) n expansion) to the standard seesaw equation for the light neutrino mass matrix, 

m v = -m T D m^m D , (4) 
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as well as the simpler relation for the heavy mass eigenstates, tun = tum- Since we are 
working in a basis where Tom is diagonal, the heavy eigenstates are then given by 

m^ ag = m N = di&g(m Nl ,m N2 ,m N3 ) . (5) 

The matrix m u can be diagonalised by the Maki-Nakagawa-Sakata unitary matrix Umns [36, 
37], leading to the following masses for the light physical states 

m^iag = U^ NS m u U MNS = diag(m I/1 ,77i I<1 ,m vs ). (6) 

Here we use the standard parameterisation for Umns given by 

( C12C13 S12C13 s 13 e~ lS \ 

Umns = s 12 c 23 - c 12 s 23 s 13 e lS c 12 c 23 - s 12 s 23 s 13 e lS s 23 c 13 V , (7) 
\ S12 S23 - C12 c 23 si 3 e J<5 -C12 s 2 3 - S12 c 23 si 3 e i<5 c 23 c i3 / 

V = diag(e-^,e-^,l), (8) 



with 



and Cij = cos 9ij, s^- = sin%. % are the neutrino flavour mixing angles, <5 is the Dirac phase 
and 0i )2 are the Majorana phases. 

In view of the above, the seesaw equation (4) can be solved for m D as [10] 



m D = t^m^R\/mt ias U MNS , (9) 

where R is a generic complex orthogonal 3x3 matrix that encodes the possible extra neutrino 
mixings (associated with the right-handed sector) in addition to the ones in Umns- can 
be parameterised in terms of three complex angles, 0j (i — 1, 2, 3) as 

/ c 2 c 3 -Ci s 3 - s 1 s 2 c 3 Si s 3 - Ci s 2 c 3 \ 
-R = c 2 s 3 cic 3 -sis 2 s 3 sic 3 - c 1 s 2 s 3 , (10) 

\ S2 S X C 2 CiC 2 / 

with q = cos^j, Si = sin^j. Eq. (9) is a convenient means of parameterising our ignorance 
of the full neutrino Yukawa couplings, while at the same time allowing to accommodate 
the experimental data. Notice that it is only valid at the right-handed neutrino scales 
Tom, so that the quantities appearing in Eq. (9) are the renormalised ones, TO^ iag (m M ) and 
Umns (^m)- 

We shall focus on the simplest scenario, where both the heavy and light neutrinos are 
hierarchical, and in particular we will assume a normal hierarchy, 

m Nl < tojv 2 < TO7V3 , 

m Ul < m V2 < m V3 . (11) 
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The masses m V23 can be written in terms of the lightest mass m Vx , and of the solar and 
atmospheric mass-squared differences as 

ml 2 = Am s 2 ol + m 2 Ui , 

ml 3 = Aml tm + m\ . (12) 

Moving now to the s-spectrum, we encounter an enlarged slepton sector due to the 
inclusion of the SUSY partners of namely v^.. The soft-SUSY breaking Lagrangian 
for the slepton sector will also include new terms and parameters. In addition to the left- 
and right-handed soft-breaking masses, mi, rug, and trilinear couplings, Ai, for the charged 
sleptons, one now includes soft-breaking sneutrino masses, m^, sneutrino trilinear couplings 
A v , and the new bilinear parameter, Bm- 

The universality conditions of the soft-SUSY breaking parameters at the high-energy 
scale M x (with M x S> tum) for the full slepton sector then read as follows: 

(m L ) I = (m^l = (mfifij = M$ % , 

( A i)ij = A (Yi)ij , = A {Y v ) i3 , (13) 

where M and A are the universal scalar soft-breaking mass and trilinear coupling of the 
CMSSM, and i,j denote lepton flavour indices, with i,j = 1,2,3. Here, we choose M x to 
be the gauge coupling unification scale. The extended CMSSM is further specified by the 
universal gaugino mass, Mi/ 2 , the ratio of the Higgs VEVs, tan (3, and the sign of the bilinear 
/^-parameter (sign//). 

The CMSSM predictions for the low-energy parameters are obtained by solving the full 
renormalisation group equations (RGEs), which must now include the appropriate equations 
and extra terms for the extended neutrino and sneutrino sectors. Due to the existence of 
intermediate scales m M introduced by the seesaw mechanism, the running must be carried 
in two steps. The full set of equations is first run down from M x to uim- At the seesaw 
scales, the right-handed neutrinos as well as their SUSY partners decouple, and the new 
RGEs (without the equations and terms for vr and vr) are then run down from uim to the 
EW scale, where the couplings and mass matrices are finally computed. 

2.2 Lepton flavour violating decays 

In the present study, and since we work within the CMSSM, all LFV originates solely from the 
neutrino Yukawa couplings. For the LFV process that we are interested in, the flavour mixing 



7 



in the neutrino sector is transmitted to the charged lepton sector via radiative corrections 
involving Y u . These corrections can be important since, due to the Majorana nature of the 
neutrinos, the Yukawa couplings may be sizable (as large as (9(1)). In particular, we will 
consider here the following LFV muon and tau decays: /i — > ej, r — > /ry, r — > ej, /i — > 3 e, 
r — > 3 e and r — > 3 /i. 

Under the requirement that at the seesaw scales {m M ) Y u satisfies Eq. (9), the running 
from Mx down to the EW scale will induce flavour mixing in the low-energy charged slepton 
squared mass matrix, M~ , whose LL, RR, LR and RL elements are given by 

M LL = + V l { Y l Yl ) t] + ^ C ° S2 P [~\ + Sin2 ^) S H> 

M RR = m l,ij + V l { Y l Yl ) l3 ~ ^ 008 2 Sin2 ^%, 

^ = (M%y , (H) 

with mz the Z-boson mass and 9w the weak mixing angle. Below rriM, the right-handed 
sneutrinos decouple, and the low-energy sneutrino mass eigenstates are dominated by the 
components [38]. Thus, the sneutrino flavour mixing is confined to the left-handed sector, 
and described by the following 3x3 matrix: 

M~ j2 = m|.. + ^m|cos2/^. (15) 

The physical masses and states are obtained by diagonalising the previous mass matrices, 
leading to 

M 2diag = R l M 2 R l\ = ^(^.^2^ 

Mf ias = W Ml ift = diag (ml, m\, m\ z ) , (16) 
where R l,v are unitary rotation matrices. 

The LFV ratios for the decay processes of our interest are obtained here via a full one-loop 
computation (and in terms of physical eigenstates), including all relevant SUSY diagrams. 

For the radiative decays lj — > lij (i ^ j), the branching ratios are given by 

BR(/,^/ l7 ) = ^-^L (|^| 2 + \A*\ 2 ) , (17) 
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where Ti. is the total lepton width, and the form factors A^ R receive contributions from two 
types of diagrams, sneutrino-chargino loops and charged slepton-neutralino loops. These 
BRs were computed in Ref. [5,6]. We will use in our analysis the explicit formulae for the 
\A\' R \ form factors as in Ref. [18]. 

Regarding the LFV decays into three leptons, lj — > 3l i: the one-loop computation was 
presented in [6], and later revised and completed in [18]. The latter work included the full set 
of SUSY one-loop contributing diagrams, namely photon-, Z-, and Higgs-penguins, as well 
as box diagrams. As shown by the explicit computation of [18], the dominant contribution 
is clearly coming from the photon-penguin diagrams and, more specifically, from the same 
form factors A\ ,R as in the case of the radiative decays 1 . This is valid even in the case of 
very large tan (3 where the Higgs-penguin diagrams, although enhanced, induce contributions 
which are still many orders of magnitude below those associated with the photon-penguins 2 . 
Therefore, the BR for the lj — > 3 k decays can be approximated by the simple expression, 

BR(lj - 3/,) = ( log^f - ^ ) x BR(Z j - , (18) 

where a is the electromagnetic coupling constant. However, and although we have verified 
that this is indeed a very good approximation, we use in the present analysis the full one-loop 
formulae for the BR(/j — > 3k) from Ref. [18]. 

Finally, and regarding the estimation of the low-energy parameters, we consider the full 
2-loop RGE running, except for the neutrino sector which is treated at the 1-loop level. 
Nevertheless, for the forthcoming discussion, it will be clarifying and interesting to compare 
the full results with the simplified estimation which is obtained within the leading logarithmic 
approximation (LLog). In the latter framework, the RGE generated flavour mixing in the 
slepton sector is summarised by the following logarithmic contributions, 

(Am|k = ~ (3 Ml + Al) (Fj L , 
(A4) y = --^AoY^YjLY^, 

(Am%) tJ — ; L\~x = log ( 5 kl , (19) 

\ m M k J 

which are originated by the running from Mx to the right handed mass scales rriMi, i = 
1,2,3. The matrix elements (YJLY,,)^ in Eq. (19) can be simply written in terms of the 

1 This has also been concluded in a generic, non-seesaw, MSSM scenario [39]. 

2 Notice that the Higgs-penguin contribution could only be relevant in a generic MSSM framework [40, 41]. 
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parameterisation of Eqs. (5-10). In particular, we obtain 



21 



£33 m Na 



Cl3 



(e* 01 jm^c X2 s 2 + e^ 2 jm^c 2 s 1 s 12 ^j + e lS ^/m^d c 2 s i3 ] 



v /m^"cic 2 ci 3 S23 -e"^ 1 v /m^"s 2 (c 23 s i2 + e l<5 c i2 s i3 s 23 ) 
+ e~^ 2 ^/m^c 2 s 1 (ci 2 c 23 - e a s i2 s 13 s 23 ) 
+ L 2 2m N2 e^ 1 ^/m^c 2 c 12 c 13 s 3 + e^ <p2 y/m^c 13 (ci c 3 - si s 2 s 3 ) s i2 
-e l<5 v /m^" (c 3 si + ci s 2 s 3 ) S13] 

- v /m^"ci 3 (c 3 si + ci s 2 s 3 ) s 23 - e~^ 1 v /m^"c 2 s 3 (c 23 s i2 + e a c i2 s i3 s 23 ) 
+e~^ 2 ^/m^ (cic 3 - sis 2 s 3 ) (c 12 c 23 - e tS s 12 s 13 s 23 ) 
+ L u m Nl e^ 1 y/m^ c 2 c 3 c 12 c 13 - e^ 2 ^/m^ c 13 (c 3 s 1 s 2 + c x s 3 ) s i2 
+e a v/m^" (- (ci c 3 s 2 ) + si s 3 ) Si 3 ] 

^frn^c 13 (- (cic 3 s 2 ) + sis 3 ) s 23 - e~^ 1 v /m^~c 2 c 3 (c 23 s 12 + e i<5 c i2 s i3 s 23 ) 
-e~^ 2 v / m^ (c 3 si s 2 + ci s 3 ) (c i2 c 23 - e J,5 si 2 s 13 s 23 ) . (20) 

The above is the relevant matrix element for the /x — >■ and /i — > 3e decays, which will be 
the most emphasised in the present work. Correspondingly, the expression for the (Yj L Y v ) 32 
((Yj LY u ) 3 i) matrix element, omitted here for brevity, is the relevant one with respect to 
r — > /i7 and r — > 3 /x (r — > e7 and r — > 3 e) decays. 

Since the dominant contribution to the fi — > e7 decay stems from the RGE induced flavour 
mixing in (Am|) 2 i, within the framework of the mass insertion and leading logarithmic 
approximations, one then obtains a simple formula given by 

a3tai1 ^ ^(3M 2 + ^) WLY„) 



BR(/i^e 7 ) = 



m SUSY 



8tt 2 



21 



(21) 



where Gf is the Fermi constant, (Yj LY V ) 21 has been given in Eqs. (20, 10), and msusY 
represents a generic SUSY mass. 



2.3 Thermal leptogenesis and gravitino constraints 

In our analysis, we will take into account constraints on LFV from the requirement of suc- 
cessfully generating the baryon asymmetry of the Universe via thermal leptogenesis [4]. In 
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this scenario, the BAU is explained by the out-of-equilibrium decay of the same heavy right- 
handed neutrinos which are responsible for the suppression of light neutrino masses in the 
seesaw mechanism. The needed CP asymmetry for BAU is obtained from the CP violat- 
ing phases in the complex angles Oi (see Eqs. (9), (10)), which also have a clear impact 
on the LFV rates. Here we assume that the necessary population of right-handed neutri- 
nos emerges via processes in the thermal bath of the early Universe. We will furthermore 
consider cosmological constraints on the reheat temperature after inflation associated with 
thermally produced gravitinos. The reheat temperature, Trh, has a strong impact on ther- 
mal leptogenesis since the thermal production of right-handed neutrinos N\ is suppressed if 
Trh < m Nl . 

2.3.1 Gravitino problems and the reheat temperature 

Thermally produced gravitinos can lead to two generic constraints on the reheat tempera- 
ture [42]. Both are associated with the fact that in the scenarios under consideration, and 
assuming R-parity conservation, the gravitinos will ultimately decay in the lightest super- 
symmetric particle (LSP). Firstly, they can decay late, after the Big Bang nucleosynthesis 
(BBN) epoch, and potentially spoil the success of BBN. This leads to upper bounds on the 
reheat temperature which depend on the specific supersymmetric model as well as on the 
mass of the gravitino. In particular, with a heavy gravitino (roughly above 100 TeV), the 
BBN constraints can be nearly avoided. In our study, we will consider the gravitino mass 
as a free parameter, so that we can safely avoid the latter constraints for any given reheat 
temperature. Secondly, the decay of a gravitino produces one LSP, which has an impact on 
the relic density of the latter. The number of thermally produced gravitinos increases with 
the reheat temperature, and we can estimate the contribution to the dark matter (DM) relic 
density arising from non-thermally produced LSPs via gravitino decay as [42] 

0-054 (^)(^). (22, 

which depends on the LSP mass, mLSP, as well as on the reheat temperature Trh- Taking 
the bound Jl£gp _th /i 2 < Qdm'* 2 i$ 0.13 from the Wilkinson Microwave Anisotropy Probe 
(WMAP) [43], we are led to an upper bound on the reheat temperature of 

r RH < 2 .4xlO-Gev(H^XY (23) 

For the considered SUSY scenarios, the mass of the LSP (which is the lightest neutralino) 
is in the range 100 GeV — 150 GeV, resulting in an estimated upper bound on the reheat 
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temperature of approximately T RH < 2 x 10 10 GeV. In the following subsection, we will 
consider the constraints on the i?-matrix parameters from the requirement of generating 
the BAU via thermal leptogenesis, while taking into account the latter bound on the reheat 
temperature. 



2.3.2 Thermal leptogenesis 

In the chosen scenario of hierarchical right-handed neutrinos, the baryon asymmetry arises 
from the out-of-equilibrium decay of the lightest right-handed neutrino iV\. The produced 
lepton asymmetry is then partially transformed into a baryon asymmetry via sphaleron 
conversion. In the MSSM, the resulting baryon to photon ratio from thermal leptogenesis 
can be written as [44] 

— « -1.04 x l(T 2 £i?7 , (24) 

n 7 

where S\ is the decay asymmetry of N\ into Higgs and lepton doublets and r\ is an efficiency 
factor for thermal leptogenesis, which can be estimated by solving the Boltzmann equations. 
The efficiency strongly depends on the ratio m Nl /Tjm as well as on the parameter rh\ [45], 
which is defined as 

m 1 = — . (25) 

In the following, regarding the estimation of the efficiency r)(rhi, TOjvi/Trh), we will use 
the numerical results of Ref. [44] for 10~ 7 eV < rh\ < 1 eV and 0.1 < j Trh < 100 
(under the assumption of a zero initial population of Ni). As presented in [44], the efficiency 
dramatically drops if either ^> Trh or if rh\ strongly deviates from its optimal value 
rh\ ~ 10~ 3 eV. Thus, the optimisation of this efficiency factor, to obtain a successful BAU, 
suggests that m Nl < T RH , which we will assume for the forthcoming LFV analysis. 

With respect to the decay asymmetries we will use the 1-loop results [46] 



In 



Xj \ Xj 



Xj 



(26) 



with Xj = m? N ./m 2 Nv for j ^ 1. 

Since in our analysis we use the i?-matrix parameterisation of Eq. (9), it is convenient 
to rewrite both the decay asymmetry e± (in the limit of hierarchical right-handed neutrinos) 
and the washout parameter mi, in terms of the i?-matrix parameters [47], 

6 m Nl E,<Im( J R? J ) 

87r V 2 J2i m ^\ R li\ 2 



-j - ~— — '—rmr- ■ ^i = ^^|%| 2 - (27) 
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Figure 1: Constraints on the i?-matrix angle 6 2 (in radians) defined in Eq. (10), from the requirement of a 
successful BAU via thermal leptogenesis (compatible with the constraint on Trh from Eq. (23)). From out- 
to inner-most rings, the regions are associated with the following BAU ranges: € [1CP 10 , 5.9 x 10 -10 ], 

n B /n 7 € [5.9 x KT 10 , 6.3 x 1CT 10 ], n B /n y e [6.3 x 1CT 10 , 1CT 9 ] and n B /n 7 > 1CT 9 . 



As seen from the previous equation, a successful leptogenesis requires complex values of the 
i?-matrix entries in order to generate a non-zero decay asymmetry. 

The BAU estimate in Eq. (24) should be compared with the reported WMAP 68% 
confidence range for the baryon-to-photon ratio [43] 

— = (6.0965 ± 0.2055) x 10~ 10 . (28) 

Finally, the constraints on the i?-matrix parameters from the requirement of a successful 
BAU compatible with the upper bound Trh ^$ 2 x 10 10 GeV are summarised in Figs. 1 
and 2. Figure 1 illustrates the impact of 82 (with 61 = 63 = 0) on the estimated BAU. As 
one can see, the 68% WMAP confidence range of Eq. (28) corresponds to a very narrow ring 
(represented by the darkest region in Fig. 1) in the Re(6 l 2)-Im(^2) plane. Notice also that 
values of either Re(6 , 2 ) or Im(# 2 ) larger than 1.2 radians (mod 7r) lead to very small values 
of the BAU,namely n^jn^ < 10~ 10 . On the other hand, the analogous study of Fig. 2 shows 
that with just 83 (61 = 62 = 0) one cannot accommodate the WMAP range. Similarly, 
values of Re(# 3 ) or Im(# 3 ) larger than 1.2 radians (mod it) also lead to excessively small 
ub/tl^ (< 10~ 10 ). Similar results regarding the constraints on 9 2 and 9 3 from successfull 
BAU via leptogenesis have been found in [17] and [19]. We also see from Figs. 1 and 2 that a 
significant part of the parameter space is excluded since the baryon asymmetry is produced 
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Figure 2: Constraints on the i?-matrix angle 6*3 defined in Eq. (10), from the requirement of successful BAU 
via thermal leptogenesis with the constraint Trh < 2 x 10 10 GeV. Colour code as in Fig. 1 (in this case 
only the region n^,jn 1 6 [10~ 10 , 5.9 x 1CP 10 ] is present). 

with the wrong sign, ris/n^ < 0, which contradicts observation. 

Regarding 81, and even though it cannot independently account for a successful BAU, 
it may have an impact on leptogenesis if 82 and/or #3 are non-zero, as can be inferred from 
Eq. (27). 

For the present study of LFV, we adopt a conservative approach, and we only require 
the estimated baryon-to-photon ratio to be within the range 

> G [l(T 10 , 1(T 9 1 . (29) 

This broad range for n^/n^ reflects the theoretical uncertainties in our estimate which may 
come, for instance, from flavour effects in the Boltzmann equations [48-50] and, more gen- 
erally, from the approximations made in [44] in order to calculate the efficiency factor 77. To 
accommodate the extended range of Eq. (29), Figs. 1 and 2 suggest that one should take 
values of 82 and 83 not larger than approximately 1 radian (mod ir). 

2.4 Implications for charged lepton EDMs 

The presence of CP violating phases in the neutrino Yukawa couplings has further implica- 
tions on low-energy phenomenology. In particular, RGE running will also induce, in addition 
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to the LFV decays, contributions to flavour conserving CP violating observables, as is the 
case of the charged lepton EDMs. Here, we also analyse the potential constraints on the 
SUSY seesaw parameter space arising from the present experimental bounds [51] on the 
EDMs of the electron, muon and tau. 

As argued in [24,25,52,53], the dominant contributions to the EDMs arise from the 
renormalisation of the charged lepton soft-breaking parameters. In particular, the EDMs 
are strongly sensitive to the non- degeneracy of the heavy neutrinos, and to the several CP 
violating phases of the model (in our case, the complex i?-matrix angles). In the present 
analysis we estimate the relevant contributions to the charged lepton EDMs, taking into 
account the associated one-loop diagrams (chargino-sneutrino and neutralino-slepton medi- 
ated), working in the mass eigenstate basis, and closely following the computation of [54, 55] 
and [23]. Instead of conducting a detailed survey, we only use the EDMs as a viability con- 
straint, and postpone a more complete study, including all phases, to a forthcoming work. 
The discussion of other potential CP violating effects, as for instance CP asymmetries in 
lepton decays, is also postponed to a future study. 

3 Results and discussion 

In this section we present the numerical results for the LFV branching ratios arising in 
the SUSY seesaw scenario previously described. In particular, we aim at investigating the 
dependence of the BRs on the several input parameters, namely on 6i, m^., and m ui , and 
how the results would reflect the impact of a potential #13 measurement. In all cases, we 
further discuss how the requirement of a viable BAU would affect the allowed parameter 
range, and in turn the BR predictions. 

Regarding the dependence of the BRs on the CMSSM parameters, and instead of scanning 
over the full (M1/2, Mo, Aq, tan/3, sign/x) parameter space, we study specific points, each 
exhibiting distinct characteristics from the low-energy phenomenology point of view. We 
specify these parameters by means of the "Snowmass Points and Slopes" (SPS) cases [56] 
listed in Table 1. 



These points are benchmark scenarios for an mSUGRA SUSY breaking mechanism. Points 
la and lb are "typical" mSUGRA points (with intermediate and large tan/5, respectively), 
lying on the so-called bulk of the cosmological region. The focus-point region for the relic 
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SPS 


Mi /2 (GeV) 


M (GeV) 


A (GeV) 


tan (3 


fi 


la 


250 


100 


-100 


10 


> 


lb 


400 


200 





30 


> 


2 


300 


1450 





10 


> 


3 


400 


90 





10 


> 


4 


300 


400 





50 


> 


5 


300 


150 


-1000 


5 


> 



Table 1: Values of Mi/ 2 , M , A , tan/3, and sign fi for the SPS points considered in the analysis. 

abundance is represented by SPS 2, also characterised by a fairly light gaugino spectrum. 
SPS 3 is directed towards the coannihilation region, accordingly displaying a very small 
slepton-neutralino mass difference. Finally, SPS 4 and 5 are extreme tan (3 cases, with very 
large and small values, respectively. Since the LFV rates are very sensitive to tan/9, we 
will also display the BR predictions as a function of this parameter, with My 2 , M , A Q and 
sign fx as fixed by the SPS points. 

To obtain the low-energy parameters of the model (and thus compute the relevant phys- 
ical masses and couplings), the full RGEs (including relevant terms and equations for the 
neutrinos and sneutrinos) are firstly run down from Mx to thm- At the seesaw scale 3 (in 
particular at mjv 3 ), we impose the boundary condition of Eq. (9). After the decoupling of 
the heavy neutrinos and sneutrinos, the new RGEs are then run down from m Nl to the EW 
scale, at which the observables are computed. 

The numerical implementation of the above procedure is achieved by means of the public 
Fortran code SPheno2 .2.2 [58]. The value of Mx is derived from the unification condition of 
the SU(2) and U(l) gauge couplings (systematically leading to a value of Mx very close to 
2 x 10 16 GeV throughout the numerical analysis), while is derived from the requirement 
of obtaining the correct radiative EW symmetry breaking. The code SPheno2 .2.2 has been 
adapted in order to fully incorporate the right-handed neutrino (and sneutrino) sectors, as 
well as the full lepton flavour structure [18]. The computation of the LFV branching ratios 
(for all channels) has been implemented into the code with additional subroutines [18]. 
Likewise, the code has been enlarged with two other subroutines which estimate the value 
of the BAU, and evaluate the contributions to the charged lepton EDMs. 

The input values used regarding the light neutrino masses and the £/mns matrix elements 

3 In our analysis we do not take into account the effect of the heavy neutrino thresholds [57]. We have 
verified that, within the LLog approximation, these thresholds effects are in general negligible in our analysis. 
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LFV process 


Present bound 


Future sensitivity 


BR(/i -> e7) 


1.2 x 10- 11 


1.3 x 10~ 13 


BR(r -> e 7 ) 


1.1 x KT 7 


10~ 8 


BR(r -> ^7) 


6.8 x 10~ 8 


10- 8 


BR(/i -> 3e) 


1.0 x 10~ 12 


10- 13 


BR(r -> 3e) 


2.0 x 10~ 7 


io~ 8 


BR(r -> 3/i) 


1.9 x 10~ 7 


10- 8 



Tabic 2: Present bounds and future sensitivities for the LFV processes. 

are 

Am 2 ol = 8 x 10~ 5 eV 2 , Am 2 tm = 2.5 x 10" 3 eV 2 , 

#12 = 30°, #23 = 45°, 013 < 10°, 5 = 0! =02 = 0, (30) 

which are compatible with present experimental data (see, for instance, the analysis of [59- 
61]). As previously mentioned, we do not address the impact of non- vanishing Umns phases 
(Dirac or Majorana) in the LFV branching ratios. The effects of Majorana phases on the 
BRs have been discussed in Ref. [26]. 

Regarding charged lepton EDMs, we require compatibility with the current experimental 
bounds [51] 

\d e \ < (6.9 ±7.4) x l(T 28 e.cm, \d^\ < (3.7 ±3.4) x 10~ 19 e.cm , 

-2.2 x 10~ 17 < d T <4.5 x 10" 17 e.cm. (31) 

Finally, and before beginning the numerical analysis and discussion of the results, we 
briefly summarise 4 in Table 2 the present LFV bounds [62-66] , as well as the future planned 
sensitivities [67-70], for the several channels under consideration 5 . 



4 In Table 2, the future prospects should be understood as order of magnitude conservative estimates of 
the projected sensitivities. 

5 There are other LFV processes of interest, such as r — > \l e e, t — > e fifi, semileptonic r decays, and 
/i-e conversion in heavy nuclei, which are not considered in the present work. With the advent of the 
PRISM/PRIME experiment at J-PARC [71, 72], /i-e conversion in heavy nuclei as Ti may become the most 
stringent test for muon flavour conservation. 
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3.1 Sensitivity to 6*13 in the case R = 1 



We begin our study by revisiting the R = 1 case which represents the situation where there 
are no further neutrino mixings in the Yukawa couplings other than those induced by the 
Uuns- I n this case, the BR(/i — * e / ~f) dependence on #13 was first observed in the context of 
SUSY GUTs [15]. In Ref. [18], a comprehensive study of all the leptonic decay channels was 
performed (in a full RGE approach), and it was noticed that /i->e7 and fi — > 3e were the 
channels that exhibited both a clear sensitivity to #13 and promising prospects from the point 
of view of experimental detection. Here, we complete the study of [18], also analysing the 
other LFV channels. More concretely, we investigate how sensitive to #13 the BR(7j — > U 7) 
and BR(/j — > 3 U) are. We also add some comments on the comparison between the full and 
the LLog approximation results. 

In Figs. 3 and 4 we plot the branching ratios of the decays fi ^ e'-f, \i — > 3e, t ^ 
c-f and r — > 3e, as a function of #13, which we vary 6 in the range [0°, 10°]. We also 
display, for comparison, the lines associated with the present experimental bounds and future 
sensitivities. In each case, we consider as input the six SPS points, and take Qi = 0, so that 
in this case no BAU is generated and there is no flavour mixing arising from the right- 
handed neutrino sector. Regarding the neutrino masses, we have assumed m Ul = 10~ 5 eV, 
while the masses of the heavy right-handed are set to = (10 10 , 10 11 , 10 14 ) GeV. In 
particular we have chosen rajvi to avoid the gravitino problem in relation with non-thermal 
LSP production, as explained in Section 2.3.1. Notice that our choice of m^ 3 leads to large 
values for the third family Yukawa couplings 7 , specifically (1^)33 ps (Y jy ) 32 ps 0.3. 

The first conclusion to be inferred from Fig. 3 is that, in agreement with [18], the sen- 
sitivity to #13 is clearly manifest in the \x — > e 7 and fi — > 3 e channels. In addition, Fig. 4 
shows that BR(r — > cy) and BR(r — > 3e) also display a strong dependence on # 13 . Notice 
that in these tau decays the BR predictions for the explored #13 values lie below the present 
and future experimental sensitivities 8 . 

The observed qualitative behaviour with respect to 13 can be easily understood from 
Eq. (20), which predicts that the dominant contribution proportional to (L 33 m^) 2 

6 The scan step is purposely finer for small values of #13. 

7 Other approaches, for instance in GUT-inspircd frameworks, allow to derive the values of mjv 3 from 
unification of the Yukawa couplings of the third family, and this may lead to even larger values of (1^)33. 
For example, an SO(10) GUT could lead to 10 15 GeV, as implied by (1^)33 « 1 (see, for example, [15]). 

8 On the other hand, we remark that compared to 13 , the uncertainties in the other neutrino oscillation 
parameters, #23, #12, Am^ and Am^ 2 , are expected to have only a smaller effect on the LFV ratios (see 
e.g. [73]) 
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Figure 3: BR(/i — > e7) and BR(^i — > 3e) as a function of #13 (in degrees), for SPS la (dots), lb (crosses), 
2 (asterisks), 3 (triangles), 4 (circles) and 5 (times). A dashed (dotted) horizontal line denotes the present 
experimental bound (future sensitivity). 
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Figure 4: BR(t 67) and BR(r — > 3e) as a function of $13 (in degrees) for SPS la, lb, 2, 3, 4 and 5. Line 
and colour codes as in Fig. 3. 



should grow as (C13 S13) 2 . For small values of #13 , the "dip" exhibited by the BRs is a 
consequence of a shift in # 13 arising from RGE running, changing it from 9 13 = 9 13(771 z) 
to 9i 3 (m M ). Renormalisation induces, in our example, that 9i 3 (m M ) m 9 Vi {ra z ) — 0.2°, so 
that the minimum of the BR is shifted from #13 = 0° to #13 ~ 0.2° (which is consistent 
with analytical estimates [74]). More explicitly, even when starting with a value #13 = 0° at 
the EW scale, RGE running leads to the appearance of a negative value for #13(7??^) (or, 
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equivalently, a non-zero positive #i 3 and 5 — 7r). 

Concerning the r — > /i 7 and r — > 3 /i channels, the corresponding branching ratios do 
not exhibit any noticeable dependence on 9 13 , as expected from the analytical expressions of 
the LLog approximation. For the case R = 1, and taking for example #13 = 5°, these BRs 
are presented in Table 3. 



BR 


SPS la 


SPS lb 


SPS 2 


SPS 3 


SPS 4 


SPS 5 


T -> /!7 


4.2 x KT 9 


7.9 x KT 9 


1.8 x KT 10 


2.6 x KT 10 


9.7 x KT 8 


1.9 x 10" 11 


r — > 3 /i 


9.4 x KT 12 


1.8 x 10" 11 


4.1 x 10~ 13 


5.9 x KT 13 


2.2 x KT 10 


4.3 x 10" 14 



Table 3: Predictions for the BR(r — > /i 7) and BR(r — > 3 /x) corresponding to the SPS points. The values of 
mjVi and are as specified in Figs. 3 and 4. In each case, the predicted values should be compared with 
the present bounds (future prospects) BR(r -> ^7) < 6.8 x KT 8 (KT 8 ) and BR(r -> 3//) < 1.9 x KT 7 

(10- 8 ). 



The conclusion to be inferred from Figs. 3, 4 and Table 3 is that, for the assumed value of 
m ui , and for the chosen seesaw scenario (which is specified by 9i and rajvj, the experimental 
bounds for BR(r — > ^7) already disfavour the CMSSM scenario of SPS 4 (for any value 
of #13). From the comparative analysis of the ^13-sensitive channels it is also manifest that 
/i — > e 7 and /x — > 3 e are the decays whose BRs are within the reach of present experiments, 
thus potentially allowing to constrain the values of $13. In fact, BR(/i — > e7) suggests that 
SPS 4, l(a and b), 3, 2 and 5 are disfavoured for values of #13 larger than approximately 0.5°, 
1°, 4°, 5° and 6°, respectively, while a similar analysis of BR(/i — > 3e) would exclude 13 
values above 1°, 3° and 4° for SPS 4, la and lb, correspondingly. Nevertheless, it is crucial 
to notice that, as can be seen from Eqs. (20, 21), the value of rriN 3 plays a very relevant role. 
For instance, by lowering mjv 3 from 10 14 GeV to 10 13 GeV one could have compatibility with 
the experimental bound on BR(yU — > ej) for #13 < 2° for all SPS scenarios. Moreover, in 
this case, even SPS 4 would be in agreement with the experimental bound on BR(r — > ^7). 

The relative predictions for each of the SPS points can be easily understood from the 
BRs dependence on the SUSY spectrum 9 and tan /3, which is approximately given by Eq. (21). 
However, it is worth emphasising that although the several approximations leading to Eq. (21) 
do provide a qualitative understanding of the LFV rates, they are not sufficiently accurate, 
and do fail in some regions of the CMSSM parameter space. In particular, for the SPS 5 
scenario, we have verified that the LLog predictions for the BRs arising from Eq. (19) differ 
from our results by several orders of magnitude. We will return to this discussion at a later 

9 For each SPS point, the associated spectrum can be found, for example, in [56]. 
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stage. 

As already mentioned, in the context of SUSY GUTs, the dependence of the BR(/i — > cy) 
on 9 13 for the same set of SPS points was presented in [15]. Instead of the full computations, 
the analysis was done using the LLog approximation, and the amount of slepton flavour 
violation was parameterised by means of mass insertions. In general, and even though a 
different seesaw scenario was considered, the results are in fair agreement with Fig. 3, the 
only exception occurring for SPS 5. In fact, while [15] predicts the largest BR(/i — > cy) for 
the SPS 5 case, our results of Fig. 3 show that the rates for this point are indeed the smallest 
ones. As already mentioned, this is due to the failure of the LLog for SPS 5. 

Henceforth, and in view of the fact that not only is the decay fi — > e 7 one of the most 
sensitive to #13, but it is also the most promising regarding experimental detection, we will 
mainly focus our discussion on the analysis of BR(/i — > e 7) . 

3.2 Implications of a favourable BAU scenario on the sensitivity 

to 13 

Motivated by the generation of a sufficient amount of CP asymmetry in the decay of the right- 
handed heavy neutrinos, one has to depart from the R—t case, and this will naturally affect 
the predictions for the several BRs. Nevertheless, it is worth stressing that the hierarchy 
of the SPS points regarding the relative predictions to the distinct LFV observables is not 
altered, and we also observe the same ordering as that emerging from Figs. 3, 4 and Table 3, 
namely BR 4 > BR ib > BR ia > BR 3 > BR 2 > BR 5 . 

As discussed in Section 2.3, the i?-matrix complex parameters 62 and 9% are instrumental 
in obtaining a value for the baryon asymmetry in agreement with experimental observation, 
while 81 plays a comparatively less relevant role. In what follows, we discuss how requiring 
a favourable BAU scenario would constrain the 9i ranges, and how this would reflect on the 
BRs' sensitivity to 9 13 . 

3.2.1 Influence of 2 

In view of the above, we begin by analysing the dependence of the BR(/x — > ej) on 9 2 and 
consider two particular values of #13, 9 13 = 0° ,5°. We choose SPS la, and motivated from 
the discussion regarding Fig. 1, take < \9 2 \ < 7r/4, with arg# 2 = 7r /4, 37r/8}. 

In Fig. 5, we display the numerical results, considering m Ul = 1CT 5 eV and m Ul = 
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Figure 5: BR(fi — > e7) as a function of \8 2 \, for arg(3 2 = {tt/8, 7r/4, 37r/8} (dots, times, diamonds, 
respectively) and #13 = 0°, 5° (blue/darker, green/lighter lines). We take m Ul — 10~ 5 (10 -3 ) eV, on the 
left (right) panel. In all cases black dots represent points associated with a disfavoured BAU scenario and a 
dashed (dotted) horizontal line denotes the present experimental bound (future sensitivity). 



10~ 3 eV, while for the heavy neutrino masses we take = (10 10 , 10 11 , 10 14 ) GeV. There 
are several important conclusions to be drawn from Fig. 5. Let us first discuss the case 
m vi = 10~ 5 eV. As previously mentioned, one can obtain a baryon asymmetry in the range 
10~ 10 to for a considerable region of the analysed \6 2 \ range. In particular, a deviation 
from the R = 1 case as small as for instance, t9 2 = 0.05 e 71 "/ 8 * can account for an amount of 
BAU close to the WMAP value. A wide region with larger values of |# 2 | (0.3 < |0 2 | < 0.8) can 
also accommodate a viable baryon asymmetry, as can be seen from Fig. 5. Notice also that 
there is a clear separation between the predictions of # 13 = 0° and t9 13 = 5°, with the latter 
well above the present experimental bound. At present, this would imply an experimental 
impact of (9 13 , in the sense that the BR predictions become potentially detectable for this 
non- vanishing t9 13 value. With the planned MEG sensitivity [67], both cases would be within 
experimental reach. However, this statement is strongly dependent on the assumed param- 
eters, in particular m vi . For instance, a larger value of m Vl = 10~ 3 eV, illustrated on the 
right panel of Fig. 5, leads to a very distinct situation regarding the sensitivity to # 13 . While 
for smaller values of \8 2 \ the branching ratio displays a clear sensitivity to having # 13 equal 
or different from zero (a separation larger than two orders of magnitude for |^| ^ 0.05), the 
effect of (9 13 is diluted for increasing values of |0 2 |. For \6 2 \ > 0.3 the BR(yU — ► ej) associated 
with #13 = 5° can be even smaller than for #13 = 0°. This implies that in this case, a 
potential measurement of BR(/i — > e-y) would not be sensitive to #13. 
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Moreover, m Ul also affects the BAU-favoured regions. In general, larger values of m Vl 
(still smaller than 1CT 3 eV) widen the range of \9 2 \ for which a viable BAU can be obtained. 
This can be understood from the fact that for very small (or zero) 9 2 and 9 3 (and with fixed 
mjvj, m Ul controls the size of the Yukawa couplings to the lightest right-handed neutrino, 
Ni. On the other side, these are the Yukawa couplings governing the washout parameter 
rh\ for thermal leptogenesis, as introduced in Eq. (25). For very small 9 2 and 9 3 , an optimal 
value rhi xs 1CT 3 eV can be reached for m Ul pa 1CT 3 eV (c.f. Eq. (27)), whereas smaller m Vl 
lead to suppressed leptogenesis in this case. For larger values of 9 2 and/or 9 3 , which can be 
still consistent with leptogenesis, m Ul becomes less important, since the other light neutrino 
masses m U2 and/or m v . A contribute to rh\ as well. In most of the following analysis, we will 
use m Ul pa 10~ 3 eV and enable a successful thermal leptogenesis by introducing a small 
i?-matrix rotation angle 9 2 . In what concerns the sensitivity to $13 via LFV, this is clearly 
a conservative choice since, as previously mentioned, lower values of m Vl (e.g. m Vl = 10~ 5 
eV) would lead to a more favourable situation. 

Whether or not a BAU-compatible SPS la scenario would be disfavoured by current 
experimental data on BR(/i — > ej) requires a careful weighting of several aspects. Even 
though Fig. 5 suggests that for this particular choice of parameters only very small values 
of 9 2 and #13 would be in agreement with current experimental data, a distinct choice of 
rriN 3 (e.g. rriN 3 = 10 13 GeV) would lead to a rescaling of the estimated BRs by a factor of 
approximately 10 -2 . Although we do not display the associated plots here, in the latter case 
nearly the entire \9 2 \ range would be in agreement with experimental data (in fact the points 
which are below the present MEGA bound on Fig. 5 would then lie below the projected 
MEG sensitivity). 

Regarding the other SPS points, which are not shown here, we find BRs for SPS lb 
comparable to those of SPS la. Smaller ratios are associated with SPS 2, 3 and 5, while 
larger (more than one order of magnitude) BRs occur for SPS 4. 

Let us now consider how the value of affects the amount of BAU, and thus indirectly 
the branching ratio associated to a given choice of 9 2 that accounts for a viable BAU scenario. 
In Fig. 6 we present the BR(/i — > ej) as a function of \9 2 \ for two distinct heavy neutrino 
spectra: m N = (5 x 10 9 , 10 11 , 10 14 ) GeV and m N = (5 x 10 10 , 10 11 , 10 14 ) GeV (values 
for m Nl respectively smaller and larger than what was previously considered). Regarding 
arg#2 ; we have chosen an example which represents a minimal deviation from the real case, 
arg#2 = 0.2, and set 9\ = 9% = 0. We consider SPS la, and again show both cases associated 
with 9±3 = 0° , 5°. From this figure, it can be seen that in the case = 5 x 10 9 GeV, 
only one BAU-favoured window is opened, for small values of 9 2 (0 < 9 2 < 7r/4). In contrast, 
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Figure 6: BR(^ — > ej) as a function of |02|, for SPS la, with arg# 2 = 0.2. 9 13 = 0° , 5° (blue/darker, 
green/lighter lines, respectively), and = 5 x 10 9 GeV, 5 x 10 10 GeV (crosses, circles, respectively). 
Dots represent points associated with a disfavoured BAU scenario for either = 5 x 10 9 or 5 x 10 10 
GeV and a dashed (dotted) horizontal line denotes the present experimental bound (future sensitivity). 



for rriN 1 = 5 x 10 10 GeV, a second window opens, corresponding to the mod it periodicity 
evidenced in Fig. 1 (also some additional points at very small \8 2 \ are allowed). The width 
of the I #2 1 interval for this second window shrinks with decreasing m^. In particular, for 
rriNi — 10 10 GeV (not displayed) this interval becomes extremely small. The latter effect can 
be understood from the interplay of 8 2 and on the relevant BAU parameters of Eq. (27). 
While rh\ is unchanged and as long as < Trh, the produced baryon asymmetry increases 
with m Nl . For a given value of m Nl , the disappearance of the second window associated 
with larger values of |0 2 | (ft < |^| ^ 3 7r/2), is due to a stronger washout, which leads to 
values of n^/n^ below the viable BAU range of Eq. (29). 

Finally, let us notice that the BAU-favoured ranges of 62 imply very distinct predictions 
for both the BRs, and the associated # 13 sensitivity. Even though the BRs arising from the 
second 82 window are significantly larger, in this case the sensitivity to #13 is considerably 
reduced, as is clearly manifest in Fig. 6. All the previous facts taken into account, we will 
often rely on the choice = 10 10 GeV and 8 2 = 0.05 e 2 * as a means of ensuring a viable 
BAU scenario via a minimal deviation from the R = t case. 
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Figure 7: BR(^ — > e7) as a function of \9i\, for arg#i = {0, — 7r/4,— 7r} (dots, times, diamonds, respec- 
tively) and (?i3 = 0° , 5° (blue/darker, green/lighter lines). BAU is enabled by the choice 02 — 0.05 e 2 * 
(f?3 =0). In all cases black dots represent points associated with a disfavoured BAU scenario and a 
dashed (dotted) horizontal line denotes the present experimental bound (future sensitivity). 



3.2.2 Influence of 6>i 

It has become clear from the previous analysis that a departure from the R = 1 case via 
non- vanishing values of # 2 can significantly affect the BR sensitivity to #13. Here we will 
show that 61 plays an equally important role on the present discussion. In Figs. 7 and 8 we 
display the BR(/i — > ej) as a function of \9\\, for different values of its argument. 

As already observed in [18], the effect of departing from the case R — 1 by varying 
61 leads to important additional contributions to the considered LFV decays. Here, we 
have only presented the case = 10 13 GeV, since for m^ 3 = 10 14 GeV the experimental 
exclusion line is already crossed for very small values of 9\ {\0\\ ~ 0.1). Opposed to the 62 
case, and as expected from the analytical estimates, there is little dependence of the BR 
on the choice of the lightest neutrino mass 10 . Considering the other SPS scenarios leads to 
analogous results, the only difference lying in a global rescaling of the BR(/x — > ej), and the 
discussion is similar to that regarding 2 . 

In the case of negative arguments, the influence of Q\ is shown in Fig. 7. Notice that in 

10 This dependence is only manifest for 6*13 « 0° and appears in terms proportional to mjv 2 , so that it is 
considerably suppressed. On the other hand, and as it occurred for 62, larger values of m Vl widen the range 
of 6>i for which a viable BAU scenario can be obtained. 
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all cases, for extremely small values of \6i\ (\6i\ < 0.1), we again recover for #13 = 5° BRs 
which are larger, and clearly distinguishable from the #13 = 0° case. In contrast, for a large 
(negative) arg 6^, the situation is reversed, and the predictions for BR(/i — > ej) associated to 
6i3 = 5° are actually smaller than for 6 13 = 0°. This becomes manifest when axgOi < — tt/2, 
a regime for which the BR starts decreasing with increasing \9i\. For real (and negative) 6\, 
(i.e. arg 6% = — it) the effect is such that for #13 = 5° two local minima of the BR are present 
(although one disfavoured by BAU), both with an associated value of the BR below the 
planned MEG sensitivity (for this specific choice of SPS point and seesaw parameters). These 
"dips" reflect a cancellation between the terms proportional to m U2 and m Vz (see Eq. (20)), 
which in fact is also present for #13 w 0°, albeit only for the second, BAU-disfavoured, 
value. It is worth pointing out that this apparent accidental cancellation for a specific 
choice of the i?-matrix parameters could correspond to the occurrence of texture zeros in 
the neutrino Yukawa couplings (motivated, for instance, from some flavour symmetry, or 
arising within specific seesaw models 11 ). Although not stable under RGE effects, these zeros 

11 At this point, we find it interesting to mention the connection between the sensitivity (or lack thereof) 
of the BR(/i — ► ej) to #13 in terms of the neutrino Yukawa couplings, as predicted within the framework 
of Sequential Dominance [75] models. For the case of "Heavy Sequential Dominance", 613 is predicted as 
#13 ~ 1(^)311/^1(^)321 + f((Y u )2i), where / is a function of the (Y„) 2 i couplings. If the predicted 6>i 3 is 
driven by the first term, and since the BR(/z — » 67) ~ (YJ)23rriN 3 L33(Y^)3i, then there is a direct connection 
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effectively translate into very small entries in the Yukawa couplings, which can account for 
the observed suppression of the BR [18] corresponding to the "dips" in Fig. 7. We would 
like to remark that, generically, the position and depth of these "dips" depend on the chosen 
values of all the seesaw parameters. 

In Fig. 8 we present a few examples of arg#i > 0. In this case, the discussion of the BRs 
and sensitivity to #13 is very similar to that conducted for small negative arguments. That 
is, for small \9±\ values, the predictions for the two #13 cases are clearly distinguishable. On 
the other hand, and irrespective of the argument (positive or negative), for sufficiently large 
|#i|, the lines corresponding to the cases 13 = 0° and 5° eventually meet, and thus for this 
choice of parameters the sensitivity of the BR to #13 is lost. 

Another relevant aspect to be inferred from Figs. 7 and 8 is how 9 1 affects the BAU 
predictions enabled by 9 2 . Unlike what occurs for 9 2 and #3, the role of 9\ in accounting for 
the observed BAU is somewhat more indirect. In particular, and as mentioned in Section 2.3, 
9 1 essentially deforms the favoured BAU areas in the 9 2 — 9 3 plane. For instance, and for 
the chosen BAU-enabling 9 2 value in Fig. 7, a real value of 9i larger than 0.9 leads to an 
estimated n B /n 7 which is no longer within the viable BAU range of Eq. (29). A distinct 
situation occurs for the cases arg#i = — 7r/4 and 7r/2, where the entire range successfully 
accounts for n B /n 7 within [10~ 10 , 10~ 9 ]. 

To conclude this subsection, let us add two further comments. Regarding the influence of 
#3 it suffices to mention that although relevant with respect to BAU (see Fig. 2), we have not 
found a significant BR(yU — > e 7) dependence on the latter parameter. This is a consequence 
of having the Yukawa couplings to the heaviest right-handed neutrino dominating, since a 
6*3 i?-matrix rotation leaves unchanged the couplings (Y u ) i3 . In this case, the sensitivity to 
9 13 is very similar to what was found for the R — t case. In the remaining analysis we will 
fix 9 3 = 0. 

Concerning the EDMs, which are clearly non- vanishing in the presence of complex 9 it we 
have numerically checked that for all the explored parameter space, the predicted values for 
the electron, muon and tau EDMs are well below the experimental bounds given in Eq.(31). 

between $13 and BR(/i — > e"f). This is, for example, what happens if we choose R = 1. In contrast, if the 
second term dominates the predicted 6*13 value, the connection of BR(/i — > e-f) to #13 is lost [76]. This 
corresponds to the "dip" of the BR(/j — > ej) in Fig. 7. 
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m N3 (GeV) 

Figure 9: BR(^ — > cy) as a function of m N3 for SPS la, with m Vl = 10~ 5 eV and m vi = eV 
(times, dots, respectively), and #13 = 0°, 5° (blue/darker, green/lighter lines). BAU is enabled by the choice 
9-2 = 0.05 e 0,2! (61 = $3 = 0). On the upper horizontal axis we display the associated value of (1^)33. A 
dashed (dotted) horizontal line denotes the present experimental bound (future sensitivity). 

3.3 Dependence on other relevant parameters: mjv 3 and tan/? 

Throughout the discussion regarding the dependence of the branching ratios on the i?-matrix 
complex angles, it has often been stressed that the leading contributions to the BRs were 
those proportional to the mass of the heaviest right-handed neutrino, m^- This is indeed 
the most relevant parameter. Here, and to briefly summarise the effect of mjv 3 , let us present 
the predictions for BR(/i — > ej) as a function of the latter mass, while keeping mjVi and 
m N2 fixed. We have checked that the BRs do not significantly depend on and mjv 2 , 
apart from one exception for m^ 2 , which we will later comment. The results for SPS la 
are displayed in Fig. 9. For completeness, we have included in the upper horizontal axis the 
associated value of (Y u ) 33 (with similar values being obtained for (Y u ) 32 )- 

We find from Fig. 9 that the full RGE result grows with mjv 3 in a very similar fashion 
to that predicted by the LLog approximation, i.e. m% 3 log 2 mjv 3 - It is clear that without 
a predictive theoretical framework for m^ 3 (e.g. GUT models) or indirect experimental 
evidence for the scale of the seesaw mechanism, there is a large uncertainty regarding the 
value of ?7ijv 3 • Within our chosen scenario of hierarchical heavy neutrinos (m^ <C m,N 2 
mjv 3 ) , assuming that the observed BAU is generated via a mechanism of thermal leptogenesis 
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(with m Nl > 10 9 GeV), and given the gauge coupling unification scale 12 (M x ~2x 10 16 
GeV), the natural choice for m N:i would lie in the range [10 10 GeV, 10 15 GeV]. It is obvious 
from Fig. 9 that such an uncertainty in rriN 3 translates into predictions for the BR ranging 
over many orders of magnitude. Hence, one can at most extract an upper bound on m,N 3 for 
the chosen set of input parameters. For instance, in Fig. 9, m Ns > 10 13 (10 14 ) GeV is not 
allowed by the present experimental bounds on the BR(fi — > cy) for 9 13 = 5° (0°). Notice 
that, although the sensitivity to #13 is clearly displayed in Fig. 9 (with more than two orders 
of magnitude separation of the #13 = 0° and 5° lines), without additional knowledge of mjv 3 
it will be very difficult to disentangle the several #13 cases. However, this argument can be 
reversed. This strong dependence on mjv 3 could indeed be used to derive hints on from 
a potential BR measurement. We will return to this type of considerations in the following 
section. 

It is also worth commenting on the local minima appearing in Fig. 9 for the lines asso- 
ciated with #i3 = 0°. As mentioned before, these "dips" are induced by the effect of the 
running of #13, shifting it from zero to a negative value. In the LLog approximation, the 
"dips" can be understood from Eq. (20) as a cancellation between the terms proportional 
to m N3 L 33 and m N2 L 2 2 in the limit #13 {tum) 0~ (with 1 — 9 3 — 0). The depth of 
the minimum is larger for smaller visible in Fig. 9. We have also checked that an 

analogous effect takes place when one investigates the dependence of BR(/i — > ej) on m^ 2 . 
It is only in this limit #13(7?^) — > 0~, and in the vicinity of the "dip", that can visibly 
affect the BRs. 

Regarding the other SPS points, with the exception of SPS 3 and 5, the results from 
the full RGE computation (not displayed here) are also in good agreement with the LLog 
approximation. The predicted BRs for SPS 3 are found to be larger than those of the 
LLog by a factor of approximately 3. This divergence is due to the fact that in the LLog 
approximation the effects of M 1 / 2 in the running of the soft-breaking parameters of Eq. (19) 
are not taken into account. Therefore, for low M and large Mi/ 2 (as is the case of SPS 3), 
there is a significant difference between the results of the full and approximate computations, 
as previously noted by [78, 79]. Moreover, this difference becomes more evident for low values 
of tan j3. 

The divergence of the two computations is more dramatic for SPS 5. This is shown in 
Fig. 10, where we compare the dependence of the BR(/i — > cy) and BR(r — » ^7) on mjv 3 , as 
given from the full computation, and in the leading log approximation. The latter approxi- 

12 The possibility of larger LFV effects arising from the existence of a higher energy scale, e.g. Mpi anc k, 
has been addressed by other authors. See, for instance [77]. 
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Figure 10: Prediction for BR(/i — > cy) and BR(r -> /17) as a function of m^ 3 for SPS 5, using the LLog 
approximation (upper blue lines) and full RGE (lower red lines). 6*13 = 0°, 5° (dots, diamonds, respectively). 
BAU is enabled by the choice 82 = 0.05 e 21 (9i = 63 = 0). A dashed (dotted) horizontal line denotes the 
present experimental bound (future sensitivity). 



mation over-estimates by more than four orders of magnitude the values of the BR(/i — > e 7) . 
The full RGE and LLog results diverge even more regarding the BR(r — > ^7), with a sep- 
aration that can be as large as five orders of magnitude. It is also manifest from Fig. 10 
that the qualitative behaviour of the full results with respect to mjv 3 is no longer given by 
m N 3 l°g 2 m N 3 - The reason for this divergence is associated to the large negative value of the 
trilinear coupling 13 , A . We considered other large negative values of A , in all cases leading 
to the same conclusion. Taking large positive Aq also leads to an important, albeit not as 
large, separation (for instance, three orders of magnitude for A = 1000 GeV). 

Finally, we briefly comment on the BR dependence on tan f3. As it is well known, the BRs 
approximately grow as tan 2 (3, and therefore this is also a relevant parameter. In Fig. 11, we 
plot a generalisation of the SPS points la and 4 (defined by M , M1/2, A Q and sign/x) with 
free tan/5, and present the sensitivity of the branching ratios to distinct values of #13. Again, 
as can be seen in Fig. 11, the sensitivity to 613 is clearly manifest, in the sense that for a given 
tan p the predictions for the BRs are very distinct for different #13 values. However, the tan (3 
dependence is so important that two tan j3 values, for instance 10 and 20, lead to predictions 
of the BR that diverge as much as those one obtains from the comparison of #13 = 3° and 
5° (for a fixed value of tan/?). This implies that unless the experimental range for tan/? 
is far more constrained than at present, we cannot conclude about the allowed/disallowed 

13 The effect of the sign of Aq in the failure of the LLog approximation has already been discussed in [78]. 
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Figure 11: BR(/x -> ej) as a function of tan/3, for toat = (10 10 , 10 u , 10 13 ) GeV, for SPS la (left) and 
SPS 4 (right). 613 = 0°,l o ,3 o ,5°, and 10° (crosses, times, asterisks, diamonds and dots, respectively). 
BAU is enabled by the choice 62 = 0.05 e 21 (Q\ = #3 = 0). A dashed (dotted) horizontal line denotes the 
experimental bound (future sensitivity) . Vertical shaded regions correspond to regions with spectra excluded 
by LEP data. 



#13 values from the present /j — > cy bounds. Just like as argued for mjv 3 , the strong BR 
dependence on tan (3 can be constructively used to further constrain tan (3 from a potential 
BR(/i — > e'y) measurement. We will address this topic in the following section. 



4 Experimental prospects: hints on SUSY and Seesaw 
parameters from measuring #13 and BRs 

In the previous section, we analysed how the several free parameters of the SUSY seesaw sce- 
nario affect the predictions for the BR(/i — > ej). We also emphasised how the sensitivity of 
the latter ratios to 613 can be altered by the uncertainty introduced from the indetermination 
of &i, tan (3 and, most of all, mjv 3 . The question we aim to address in this section is whether 
a joint measurement of the BRs and #13 can shed some light on apparently unreachable 
parameters, like m^ 3 . 

The expected improvement in the experimental sensitivity to the LFV ratios (see Ta- 
ble 2) support the possibility that a BR be measured in the future, thus providing the first 
experimental evidence for new physics, even before its discovery at the LHC The prospects 
are especially encouraging regarding /1 — > ej, where the sensitivity will improve by at least 
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Figure 12: Contours of BR(/i — > ej)= (1.2 ± 0.1) x 10~", with n = 10, . . . , 15 in the m,N 3 — tan/3 plane, 
for a generalised SPS la scenario. We assume 6*13 = 1° ± 0.1° (5° ± 0.5°) on the left (right) panel. BAU is 
enabled by the choice 82 — 0.05 e 21 (6>i =63 = 0). The current experimental bound is associated with the 
darkest (blue) surface, while the future sensitivity is represented by the lightest (cyan) one. 



two orders of magnitude. Moreover, and given the impressive effort on experimental neutrino 
physics, a measurement of 6^3 will likely also occur in the future [27-35]. In what follows, 
let us envisage a future "toy"-like scenario, where we will assume the following hypothesis: 
(i) measurement of BR(/i — > ej); (ii) measurement of 613 ; (iii) discovery of SUSY at the 
LHC, with a given spectrum. Furthermore, we assume that BAU is explained via thermal 
leptogenesis, with a hierarchical heavy-neutrino spectrum. 

Under the above conditions, let us conduct the following exercise. First, choosing SPS la, 
m Nl = 10 10 GeV, m N2 = 10 11 GeV, m Ul = 1(T 3 eV, 6 2 = 0.05 e°- 2i (a minimal BAU-enabling 
deviation from the R = 1 case), and with #13 set to 1° (±0.1°) and to 5° (±0.5°), we predict 
the BRs as a function of tan/? and rriN 3 - We then plot the contour lines for constant BR 
values in the — tan (3 plane. In Fig. 12 we display the corresponding contours for the 
central values of 1.2 x 10~ n with n — 10, . . . , 15, allowing for a 10% spread-out around these 
values. The predicted contours should be compared with the present bound and future 
sensitivity of 1.2 x 10 -11 [62] and 1.3 x 10~ 13 [67], respectively. 

Given a potential SUSY discovery, the implications of a measurement of BR(/i — > cy) 
and #13 are clearly manifest in Fig. 12. From this figure we first learn that, even in the 
absence of an experimental determination of tan f3, a potential measurement of BR(/i — > e 7) 
and #13 will allow to constrain m^ 3 . For example, an hypothetical measurement of BR(yU — > 
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Figure 13: Contours of BR(/i — > ej)= (1.2 ± 0.1) x 10 n , with n = 10, . . . , 15 in the tojv 3 — tan/3 plane for 
a generalised SPS 2 (left) and SPS 4 (right) scenarios. We assume 613 = 5° ± 0.5°. BAU is enabled by the 
choice f5 2 = 0.05 e a2i (6>i = 3 = 0). 

e7)~ 1.2(±0.1) x 10 -12 would point towards the following allowed ranges of m^. A : 

013 w 1° 2 x 10 13 GeV < m N3 < 2 x 10 14 GeV , 

9 n w 5° 1.5 x 10 12 GeV < < 10 13 GeV. (32) 

Other assumptions for the BRs would equally lead to an order of magnitude interval for the 
constrained values of mjv 3 . If in addition to the s-spectrum, we assume that tan/3 is exper- 
imentally determined, then the intervals for mjy 3 presented in Eq. (32) can be significantly 
reduced. For instance, assuming that SPS la is indeed reconstructed (that is, tan/3 = 10), 
then we would find 

13 « 1° 4 x 10 13 GeV < m N:i < 7 x 10 13 GeV , 

6 n w 5° 3 x 10 12 GeV < m Na < 5 x 10 12 GeV . (33) 

The hypothetical reconstruction of any other SPS-like scenario would lead to similar one 
order of magnitude intervals for mjv 3 but with distinct central values. As expected, 
for the same BR and #13 measurements, SPS 2, 3 and 5 lead to larger values of mjy 3 , the 
contrary occurring for SPS lb and 4. This can be seen in Fig. 13, where we present an 
analogous study to that of Fig. 12, but focusing on SPS 2 and SPS 4, and only considering 
#13 = 5° ±0.5. 

Concerning the comparison with current experimental bounds, one can also draw some 
conclusions regarding the excluded regions of the 777,^3 -tan /? plane. From both Figs. 12 
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and 13, for 613 = 5°, and for the chosen set of input parameters, we infer that in all cases 
the upper-right regions of the m^-tan/? plane are clearly disfavoured. For instance, for 
SPS la, m N . A > 10 14 GeV would be excluded for any value of tan/?. In the case of SPS 
2, the exclusion region would be delimited by tan/5 > 10 and m N[i > 10 14 GeV. The most 
pronounced exclusion region is for SPS 4, and is given by tan/3 > 20, m Ns > 10 13 GeV. 
With the expected future sensitivity, these exclusion regions will be significantly enlarged. 

A potential caveat to the previous discussion is the fact that, as seen in Section 3, there 
is a very important dependence of the BRs on the i?-matrix parameters 0j. Not only will 
this have implications on how accurate the indirect estimates of m Nz are, but will also affect 
any judgement regarding the experimental viability of a SUSY seesaw scenario. We recall 
that, as shown in Section 3, other choices of 9 2 (and 0-y) can lead to substantially smaller or 
larger BRs, therefore modifying the exclusion regions of Figs. 12 and 13. 

To take into account the strong i?-matrix dependence, let us conduct in what follows a 
more comprehensive survey of the parameter space. For SPS la, and for distinct choices of 
the heaviest neutrino mass, we scan over the BAU-enabling i?-matrix angles (setting 6> 3 to 
zero) as 



Given that, as previously emphasised, /x — > cy is very sensitive to 9\s, whereas this is not 
the case for BR(r — > /J7), and that both BRs display the same approximate behaviour with 
m^ 3 and tan/5, we now propose to study the correlation between these two observables. 
This optimises the impact of a #13 measurement, since it allows to minimise the uncertainty 
introduced from not knowing tan f3 and m^ 3 , and at the same time offers a better illustration 
of the uncertainty associated with the i?-matrix angles. In this case, the correlation of the 
BRs with respect to mjv 3 means that, for a fixed set of parameters, varying rriN 3 implies that 
the predicted point (BR(r — > /J 7), BR(/i — > ej)) moves along a line with approximately 
constant slope in the BR(r — > /i7)-BR(/x — ■> ej) plane. On the other hand, varying 9 13 
leads to a displacement of the point along the vertical axis. In Fig. 14, we illustrate this 
correlation for SPS la, and for the previously selected rriN 3 and #i )2 ranges (c.f. Eq. (34)). We 
consider the following values, 9 13 = 1°, 3°, 5° and 10°, and only include the BR predictions 
allowing for a favourable BAU. In addition, and as done throughout our analysis, we have 
verified that all the points in this figure lead to charged lepton EDM predictions which are 
compatible with present experimental bounds. More specifically, we have obtained values 



< < tt/4, 
< |02| < tt/4, 
m Nz = 10 12 , 10 1 
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-tt/4 < arg^ < tt/4, 
< arg# 2 < tt/4, 
10 14 GeV. 



(34) 
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Figure 14: Correlation between BR(/i —> ey) and BR(r — > ^7) as a function of mjv 3) for SPS la. The 
areas displayed represent the scan over 9i as given in Eq. (34). From bottom to top, the coloured regions 
correspond to #13 = 1°, 3°, 5° and 10° (red, green, blue and pink, respectively). Horizontal and vertical 
dashed (dotted) lines denote the experimental bounds (future sensitivities). 

for the EDMs lying in the following ranges (in units of e.cm): 
10~ 39 < 1 41 < 2xl(T 35 , 6xlCr 37 < |^| < 1.5xl0- 32 , 1(T 34 < \d T \ < 4x 1(T 31 . (35) 

For a fixed value of mjv 3 , and for a given value of #13, the dispersion arising from a 9\ 
and 62 variation produces a small area rather than a point in the BR(r — > /i7)-BR(/x — > ej) 
plane. The dispersion along the BR(r — > fi 7) axis is of approximately one order of magnitude 
for all 13 . In contrast, the dispersion along the BR(/i — > ej) axis increases with decreasing 
9i3 (in agreement with the findings of Section 3), ranging from an order of magnitude for 
9\3 = 10°, to over three orders of magnitude for the case of small 6*13 (1°). From Fig. 14 we can 
also infer that other choices of mjv 3 (for # 13 £ [1°, 10°]) would lead to BR predictions which 
would roughly lie within the diagonal lines depicted in the plot. Comparing these predictions 
for the shaded areas along the expected diagonal "corridor" , with the allowed experimental 
region, allows to conclude about the impact of a 9x3 measurement on the allowed/excluded 
mjv 3 values. 

The most important conclusion from Fig. 14 is that for SPS la, and for the parameter 
space defined in Eq. (34), an hypothetical 6^3 measurement larger than 1°, together with 
the present experimental bound on the BR(/i — ► e^), will have the impact of excluding 
values of > 10 14 GeV. This lends support to the hints already drawn from Fig. 12. 
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Figure 15: Correlation between BR(/i — > ej) and BR(r — * /17) as a function oim,N 3 - The points displayed 
mimic the behaviour of the central #13 = 10° marked point in Fig. 14, for SPS points la, lb, 2, 3, 4 and 5. 
Horizontal and vertical dashed (dotted) lines denote the experimental bounds (future sensitivities). 



Moreover, with the planned MEG sensitivity, the same #13 measurement can further constrain 
rriN 3 < 3x 10 12 GeV. The impact of any other #13 measurement can be analogously extracted 
from Fig. 14. 

Similar conclusions can be reached for the other SPS points, as seen in Fig. 15, where 
we only display the predictions corresponding to the point marked with a cross in the centre 
of the #i3 = 10° shaded area of Fig. 14 (taking into account all #13 values would lead to 
replications of the shaded areas observed in Fig. 14). Regarding SPS lb, the discussion is 
very similar to that of SPS la, and the inferred constrains on m^ 3 are almost identical. 
SPS 2 and SPS 3 offer very close predictions, and when compared to SPS la, for the same 
#13 measurement, allow to extract slightly weaker bounds on m^ A . On the other hand, 
SPS 4 clearly provides the most stringent scenario and a measurement of #13 = 10° is only 
compatible with m^ 3 < 10 12 GeV. Notice also that this is the only case where the present 
experimental bound from BR(r — > /X7) plays a relevant role. SPS 5 provides the weakest 
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bounds on m,N 3 but nevertheless still allows to exclude mjv 3 > 10 14 GeV from a measurement 
of #13 = 10°. Finally, it is interesting to notice that the observed correlations for SPS 5 are 
manifestly different from the other cases, in agreement with the findings of Section 3. In 
this case, varying mN 3 leads to predictions of the BR(/i — > ej) and BR(r — > fxj) which are 
not linearly correlated, opposed to what would be expected from the LLog approximation. 

5 Conclusions 

In this work, we have investigated lepton flavour violating muon and tau decays in the 
CMSSM extended by three right-handed (s)neutrinos, and used a type-I seesaw mechanism 
to explain the smallness of the neutrino masses. As typical examples of an mSUGRA-like 
scenario, several SPS SUSY benchmark points were considered. We have parameterised the 
solutions to the seesaw equation in terms of a complex orthogonal matrix R and of the right- 
handed neutrino masses, requiring compatibility with low-energy data. We have considered 
scenarios of hierarchical light and heavy neutrinos. In addition, we imposed consistency with 
present bounds on charged lepton EDMs and baryogenesis via thermal leptogenesis taking 
into account constraints on the reheating temperature from non-thermal LSP production by 
gravitino decay. We have studied in great detail the sensitivity of the BRs to 9 n , giving 
special emphasis to the \i — > e 7 decay channel. 

In a first stage, we have considered the simple case R = 1 where there are no additional 
neutrino mixings other than those in the Umns- We have found a very pronounced sensitivity 
to #13 in the decay channels \i — > ej, /i — *3e, r ^ e7 and r — > 3e. Varying #13 from 0° 
to 10°, the branching ratios for the above processes increase by several orders of magnitude. 
In view of the present experimental bounds and the expected future sensitivity, /i — > ej 
is by far the most promising channel to study the sensitivity to 13 in LFV processes. We 
would like to notice that /x — * 3 e may also offer interesting additional information. We have 
presented the predictions for the branching ratios for various SPS SUSY benchmark points. 
We further emphasised the importance of a full numerical computation, which we have found 
to differ significantly from the LLog approximation in some cases. 

We have then explored how the sensitivity of BR(/i — > ej) to 9 13 is altered when we take 
into account the remaining SUSY seesaw parameters. In this sense, we have found that the 
most relevant parameters are #1, 62, rriN z and tan/3 and we have systematically studied their 
influence on the impact of 9 13 on BR(/i — > e^). We have also noticed that the sensitivity to 
#13 improves for lower values of m Vl (m Ul < 10~ 3 eV). 
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Compared to the special case R—t, non- vanishing 9i can have important consequences. 
In particular, the sensitivity to 9 13 is considerably reduced for large values of \6 X \ and |6> 2 |. 
However, thermal leptogenesis severely constrains 9 2 and 9 3 (but not generically 9\). In 
fact, the requirement of successful thermal leptogenesis with constraints on the reheating 
temperature from non-thermal LSP production by gravitino decay, suggests a region |6 , 2i3 | < 
1. In these ranges, we have studied the sensitivity of the BR to 9 13 . Generically, the 
separation between the BR predictions for distinct #13 is reduced when we move from R—t 
to R 7^ 1, and one could be led to the conclusion that the BR sensitivity to 9 13 would be 
reduced. However, we have also found cases of R ^ t where, although this separation is 
reduced, the BR predictions are now larger (and can be above the experimental bounds) and 
different 6> 13 values can be distinguished even more efficiently than in the R — t case. 

Regarding the right-handed neutrino masses, the most relevant one for the LFV BRs is 
clearly m^ z (with a marginal role being played by m N2 ). Even though m Nl does not directly 
affect the BRs, it nevertheless plays a relevant role with respect to baryogenesis. This, 
together with the assumption of hierarchical right-handed neutrinos, leads furthermore to 
an indirect lower bound for rriN 3 - For a given choice of #13, the dependence on mjv 3 is so 
pronounced that for the investigated range [10 11 GeV, 10 15 GeV], the BRs change by over 
six orders of magnitude. Thus, and even though the sensitivity to 9 13 is clearly manifest 
(for instance, more than two orders of magnitude separation between the BR predictions of 
#13 = 1° and 5°, for a given value of 6> 2 ) without additional knowledge of mjv 3 it will be very 
difficult to disentangle the several 9 n cases. 

In a similar fashion, the sensitivity of the BRs to 9± 3 can be altered by the uncertainty 
introduced from the indetermination of tan f3. The study of the generalised SPS points shows 
that changing tan f3 from 5 to 50 translates in BR(/z — > e 7) predictions which differ by two 
orders of magnitude, so that unless there is an experimental determination of tan j3, it will 
also be hard to distinguish the distinct #13 predictions. Moreover, we have emphasised that 
this strong dependence on mjv 3 and tan (3 can be constructively used as a means of extracting 
information on these parameters from a potential joint measurement of #13 and BR(yU — > e 7). 

Remarkably, within a particular SUSY scenario and scanning over specific 9\ and #2 
BAU-enabling ranges for various values of #13, the comparison of the theoretical predictions 
for BR(/i — > cy) and BR(r — > ^7) with the present experimental bounds allows to set 
0i3-dependent upper bounds on mjv 3 . Together with the indirect lower bound arising from 
leptogenesis considerations, this clearly provides interesting hints on the value of the seesaw 
parameter m N3 . For instance, in the SUSY scenario SPSla and for values of 9 13 in the present 
experimental allowed range, the present MEGA constraint on BR(yU — > ej) already sets an 
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upper bound on m Ns , m Nz < 10 13 GeV for 13 > 10° and m Na < 6 x 10 13 GeV for 13 > 3°, 
as inferred from Fig. 14. These bounds are even more stringent for the case of SPS4 (see 
Fig. 15) where the present contraint on BR(/i — > cy) sets an upper bound of m Na <3x 10 12 
GeV for 13 > 10° and m Ns < 10 13 GeV for 13 > 3°. With the planned future sensitivities, 
these bounds would further improve by approximately one order of magnitude. 

Ultimately, a joint measurement of the LFV branching ratios, 6> 13 and the sparticle spec- 
trum would be a powerful tool for shedding some light on otherwise unreachable SUSY 
seesaw parameters. 
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